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An improved acceleration scheme to produce protons with controlled divergence and con¬ 
centrated energy density is studied using ultrashort ultraintense (USUI) laser pulse interac¬ 
tion with a tailored hole-target in target normal sheath acceleration (TNSA) regime. Two- 
dimension-in-space and three-dimension-in-velocity (2D3V) particle-in-cell (PIC) simulations 
show that the tailored hole-target helps to reshape the sheath electric field and generate a 
transverse quasistatic electric field of TVjm along the inner wall of the hole. The transverse 
electric field suppresses the transverse expansion of the proton beam effectively, as it tends 
to force the produced protons to focus inwards to the central axis, resulting in controlled 
divergence and concentrated energy density compared with that of a single plain target. The 
dependence of proton beam divergence and energy feature on depth of the hole is investi¬ 
gated in detail. A rough estimation of the hole depth ranges depending on uq of the incident 
laser is also given. 
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I. INTRODUCTION 

The laser-driven ion acceleration from ultrashort ultraintense (USUI) laser pulse interaction with 
various solid targets has been studied actively for applications ranging from basic particle physics 
[1], bench-top particle accelerators[2]-[4], medical therapies[3], fast ignition of inertial controlled 
fusion (ICF)[4][5], etc. Up to now, several mechanisms for accelerating ions have been proposed, 
such as target normal sheath acceleration(TNSA)[2][8], laser breakout after-burner (BOA) [9], and 
radiation pressure acceleration(PRA)[10][ll]. Many potential applications require proton and ion 
beams with high collimation, monoenergetic, larger particle number and intense energy density, as a 
result, the enhancement of beam quality becomes of intriguing interest and numerous experimental 
and theoretical studies have been devoted to achieve this goal[3][12]-[17]. 

In TNSA regime, which is more stable than other ions acceleration mechanisms, when an intense 
laser pulse irradiates on a thin plain target, a large number of electrons are accelerated by the 
incident laser and then transport to the target backside, forming an electron cloud and thus a 
strong electrostatic charge-separation field there. A population of protons near the rear surface of 
the target then are pulled out and accelerated into the backside vacuum by this so-called target 
normal sheath electric field[2][18]. One of the challenges here is how to generate a nearly local and 
uniform sheath field at the target rear surface, as typically the sheath field occupies a large area and 
is inhomogeneous due to the divergent electrons which establish the sheath field, and more to the 
point, when the normally incident laser is p-polarized, in which case the electric field of the laser 
oscillates along the transverse direction and thus pushes electrons to move up or down to either 
lateral edges of the plain target and assemble there, as a consequence, the target normal sheath 
electric field is stronger on both the lateral edges symmetrically while relatively much weaker on 
the central axis, and this edge effect will get greater over time. In order to control the shape of the 
sheath electric field and then to obtain a collimated proton beam with high quality, several tailored 
structural targets and ions-doped foil targets are proposed previously[19] [20]. 

In this paper, we study a practical scheme to generate proton beams with controlled diver¬ 
gence and concentrated energy in TNSA regime by a USUI laser illuminating the tailored hole- 
target, and we also employ a same single plain target without the backside hole as a comparison. 
Two-dimension-in-space and three-dimension-in-velocity (2D3V) particle-in-cell (PIC) simulations 
demonstrate the effectiveness of the tailored hole-target in suppressing the transverse proton beam 
divergence by confining the sheath electric filed in the hole almost locally and uniformly and also 
by generating a transverse electric field to focus the protons. Accordingly, the dependence of the 
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proton beam characteristics on the target hole depth is also worthy to investigate in detail. 

This work is organized as follows. In Sec. II, we present the target model and the simulation 
parameters. For comparison, both of the single plain and tailored hole-target are considered. In Sec. 
Ill, our simulation results are presented, from which one can see the robust improvement on the 
beam divergence and protons energy intensity by using the tailored hole-target. The dependence 
of proton beam divergence on depth of the hole is examined thoroughly in Sec. IV, as well as some 
estimations about hole depth range on laser intensity ao- A conclusion is given in Sec. V. 


II. TARGET CONFIGURATION AND SIMULATION PARAMETERS 


The simulations are performed with a 2D3V PIC code KLAP2D[17]. In the simulations, 1500 
cells longitudinally along 2 ;—axis and 2000 cells transversely along axis constitute a 15/im x 20/im 
simulation box. A p—polarized laser pulse with ao = 6 (intensity Iq ~ 4.9 x lO^^W/cm?) and 
wavelength Aq = 1/im is normally incident on a solid plain target comprised of electrons and 
protons. The front side of the plain target is located at z = 5/im, of which thickness is 1/im and 
width is 20/im. In order to include the prepulse effect, we employ a linear density gradient in 0.5/im 
at the laser illumination surface. Figure shows a conceptual diagram of the single plain target and 
the tailored hole-target studied in our simulation. In the case of the single plain target, see Figs. 
1(a), the initial electrons and protons peak density is 507ic, and the tailored hole-target is designed to 
have the same configuration and particle settings with the single plain one, and together with these, 
a hole with depth of h and diameter of 4/im surrounding by two 5/im thick, lOOric dense horizontal 
ramparts made up of and electrons, see Figs. 1(b). The initial temperature of electrons is set 
to be lOOOeU. About (1.4 ^ 2.0 x 10^) superparticles are employed in our simulations. The laser 
pulse coming from the left boundary has a transverse Gaussian profile with beam waist tq = 2.biim 
and a trapezoidal temporal profile of duration r = 22T, consisting of a plateau of 20T and rising 
and falling periods of IT each, where T is the laser period. As has been studied by T. P. Yu et 
al.|13j and M. Nakamura et al.[T9], the optimal hole diameter is of the order of the laser spot size 
and proton beam divergence and energy characteristics have little dependence on pulse duration, 
so we fix the laser parameters and hole diameter all through this paper. 
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Figure 1. The conceptual diagram of the single plain target and the tailored hole-target studied in our simulation, 
(a) the single plain target with the initial electrons and protons peak density of 50nc, and (b) the tailored hole-target, 
which is designed to have the same configutation and particle settings with the single plain one, and together with 
these, a hole with depth oi h — and diameter of 4/im surrounding by two thick, lOOric dense horizontal 
ramparts made up of and electrons. The front side of the plain target is located at z = 5/xm, and its thickness 
and width is Ifim and 20/i?Ti, respectively. In order to include the prepulse effect, we employ a linear density gradient 
in 0.5fim at the laser illumination surface.The initial temperature of electrons is set to be lOOOeF. The target 
materials have been marked. 

III. EFFECT OF HOLE-TARGET IN RESHAPING SHEATH FIELD AND INDUCING 

TRANSVERSE ELECTRIC FIELD 

The main acceleration mechanism we consider here is TNSA, which means the shape of sheath 
electric field determined by the accelerated hot electron cloud spread is vital to the quaity of the 
subsequently generated proton beam. Firstly, the shape of sheath electric fileds for both plain and 
tailored hole-target are investigated. Figs. |^a)-(f) show distributions of cycle-averaged sheath field 
Esheath {Esheath IS actually the longitudinal electric field E^) at t = lOT, 20T, SOT for plain (top) 
and hole-target (bottom), respectively. One can notice that in Figs. |^a)-(c) for the plain target, 
at early time, say t = lOT, the maximum sheath field ^ 1.2 x lO^^V/m is centralized 

around the laser incident axis with a diameter of about 10.3A and the further away from the 
central axis, the weaker Esheath- However, as the accelerated hot electrons are further exploding 
into the vacuum, say at t = 20T, Egheath is expanding along both the transverse and longitudinal 
directions as well, resulting in a wider and longer bell shape. What’s more, Esheath is still y- 
axial symmetric but is stronger ^ 6.6 x lO^^U/m) at the target top and bottom edges 

than that ^ 4.5 X lO^^U/m) in the region nearer around laser axis. At later time points 

T > SOT, Esheath has expanded more decentralized, and the maximum sheath field is obviously 
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located at the top and bottom edges of the plain target, leaving the centraxonial E sheath much more 
weaker. In addition, because the incident laser is p-polarized, the electric field of the laser pushes 
electrons along the transverse direction, and thus makes the hot electron denser at the edges. Add 
it all up, the transverse edge effect of Esheath is one of the reasons which lead to proton beam 
divergence and large spot size. As shown in Figs, [^d)-(f), in the case of the hole-target, the 
transverse edge effect of Esheath is suppressed by the horizontal ramparts made of Al, which forces 
Esheath fo distribute uniformly and locally in the hole, and from our simulation results, we obtain 
Esheath ^ 9.0 X lO^^y/m, 4.8 X lO^^y/m, 2.1 x at t = lOT, 20T, SOT, respectively. The 

reason why the hole can help to eliminate the edge effect is the accelerated hot electrons can not 
transport freely in the closely-behind Al ramparts, in which denser electrons prohibit them from 
going forward and thus Esheath can not be established in these regions like that of the plain target. 



z (g m) z (g m) z (g m) 

Figure 2. Distributions of cycle-averaged sheath electric fields Esheath at t = lOT, 20T, SOT for plain (a)-(c) and 
hole-target (d)-(f), respectively. The dashed black lines show the initial inner boundaries of the Al ramparts. The 
electric fields are normalized by ^ where rUe , cuo , c and e are electron rest mass, laser angular frequency, light 

speed in vacuum and electron charge, respectively. 

As the shape of Esheath evolves over time, we pick up three snap shots of sectional view (the 
crosssection is at z = 6.5/im ) at the same time points mentioned above, and the results are 
presented in Figure]^ Compared with Esheath of the plain target (red line), which is much broader 
and has two sharp corners almost symmetrical about the central axis (y=10/im), Esheath of the 
tailored hole-target (black line) is limited just tightly around central axis, being local and uniform. 
Moreover, for the single plain target, as we discussed above, the two corners of Esheath moving 
outwards to the top and bottom edges of the plain target and meantime the central part of Esheath 


































6 


is getting weaker. The maximum Egheath shown in Figure is smaller for the tailored hole-target 
due to less accelerated hot electrons produced in the hole than that in vacuum of the plain target 
rear side. 



Figure 3. Sectional view for the crosssection located at z = 6.5/xm at time points (a) t = lOT, (b) t = 20T, (c) 
t = SOT, red line is for the plain target and black line is for the tailored hole-target. The electric fields E sheath are 
cycle-averaged and normalized by . 


The tailored hole-target not only has an advantage over the plain one in controlling the shape 
of Esheath 1 but also has a transverse electric field Ey induced in the hole which does not exist in 
the case of the plain target. As one can see in Figure Ey has opposite directions in the hole, 
that is, in the upper side Ey is negatively along —y while in the lower side Ey is positively along 
+y^ thus this transverse field tends to focus and confine the protons in the transverse direction as 
a tight bunch. The transverse electric field can reach as large as 6.43 x lO^^y/m according to our 
simulation results, which is approximately consistent with 


Ey ~ 


eXu 


(^1 + al - l)meC^ 


eXo 


= 6.71 X lO^^V/m 


where we estimate Xd = A/eoTg/nge^ as the Debye length of accelerated hot electrons with temper¬ 
ature Tg PS rrieC^xIl + ^ 2.23MeV. 

y Thig, c 
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Figure 4. (a) Distribution of transverse electric field Ey in the hole at t = 20T; (b) sectional view of transverse 
electric field Ey for the crosssection located at z = llfim at t = 20T, where the black dashed line indicates the 
crosssection. The electric fields Ey are cycle-averaged and normalized by 
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In order to have a clear understanding about the energy distribution characteristics of generated 
protons in the rear side of the targets, we use Figure to show the distribution of proton energy 
density in the region beyond 2 ; = 5/im at time t = SOT for both plain and tailored hole-target. 
One can ensure that thanks to both the reshaped Egheath and the induced transverse electric field 
Ty, protons produced by the tailored hole-target have been confined in the hole with energy more 
greatly concentrated compared with that of the plain target. 



Figure 5. Distributions of proton energy density in the region z > 5/im at t = SOT for (a) the plain target, (b) 
tailored hole-target. Protons energy is in unit of meC^. 

IV. THE DEPENDENCE OF PROTON BEAM DIVERGENCE ON DEPTH OF THE 
HOLE AND DEPTH RANGE ESTIMATION 

In Section II when declaring the simulation parameters, we fix diameter of the hole (D = 
4/im) and leave depth h as the only geometric variable for good reason. Here we investigate the 
dependence of proton beam divergence and energy feature on h. Figure demonstrates the rear 
sheath field Egheath at t = SOT for the single plain target and tailored hole-target with different 
depth h = 0.6/im, l/im,5/xm, which are given respectively by lines in different colors as stated in 
legend. As we have proved above, there is severe edge effect for the plain target, leading to proton 
divergence. All the hole-targets have their sheath electric fields limited in the inner hole, except 
for the difference that the field distributions are distinct marginally (at hole inner margin 2 ; = 8/im 
and 2 ; = 12/im) and in the hole. In the case of /i = 0.6/im, marginal sheath field is the strongest 
among the four targets, which may result in more accelerated electrons in marginal region and thus 
large beam divergence; and further more, ’’marginal effect” makes Esheath fluctuate more wildly 
in the hole. For h = 1/im case, we still see the ’’marginal effect”, but a little eased. The most 
optimal shape of sheath electric field for accelerating protons is in the case of /i = 5/im, where 
the marginal fields are almost suppressed to equal the level of the electric field in the hole, and 






therefore, the sheath field Egheath is reshaped to be uniform and centralized. From our simulation 
results, we ensure that h = 5/im is the optimal depth parameter among these cases, for too shallow 
the ’’marginal effect” affects Esheath significantly. 



Figure 6. Distribution of rear sheath electric fields Esheath for the single plain target and tailored hole-target with 
different depth h = O.G/xm, l/xm, 5/im at t = SOT, where the black line is for the plain target, red, green and blue 
lines are for h = 0.6, /i = 1 and /i = 5, respectively. The electric fields Esheath are cycle-averaged and normalized by 

meUJQC 

e 


Now we focus on the dependence of proton beam divergence on depth of the hole. Figure 
shows the divergence spectrum of accelerated forward-going protons in the rear side at t = SOT for 
the single plain target and tailored hole-target with different depth h — 0.6/im, 1/im, 5/im. Here 
divergence angle is defined by the following formula 

^div — arc tan —^ 

Pz 

where py and pz are protons transverse and longitudinal relativistic momentum respectively. One 
can see clearly that the proton beam from the plain target has two divergence angle peaks in 
^div = 3.946° and Odiv = 1.91°, while proton beam from the hole-target with h = 0.6/im has an 
obvious angular deviation from the central axis and one angle peak in O^iv ~ 1.32°, and proton 
beams from the hole-targets with h — 1/im and h — 5/im each has angle peaks in Odiv ^ 0.04° and 
^div ~ 0.19° respectively. Although the hole-target with h = 1/im has a smaller peak divergence 
angle, on the whole its divergence spectrum shape is fatter than hole-target with h — 5/im. The 
accelerated proton number in the case h — 5pm from our simulation results is about 9.65 x 10^. 

So far, we have learned the primary function of the ramparts of the hole-target is to eliminate the 
electric sheath field ’’edge effect”, confine the accelerated hot electrons and focus the subsequent 










9 



-20 -10 0 10 20 
Divergence (degree) 

Figure 7. Divergence spectrum of accelerated forward-going protons in the rear side at t = SOT for the single plain 
target and tailored hole-target with different depth, where the black line is for the plain target, red, green and blue 
lines are for h = 0.6/i?n, h = l/xm, and h = respectively. The divergence angle is in unit of degree, and the 

proton number is indicated in the vertical axis. 


4 6 


10 


Figure 8. A rough estimation of the hole depth ranges depending on ao of the incident laser. 

protons, and we have demonstrated the dependence of proton beam divergence on depth of the 
hole, and next, the depth range should be figured out. As we know, the accelerated hot electrons 
temperature is Tg ~ meC^(l + and according to the following formula 

Tg = eEiLn 

we can estimate the minimum depth of the hole should be the local plasma scale length = Cgt ^ 
0.47(1 + where Cg = VTii is the ion sound speed, Ei is the longitudinal accelerating 

field and t is taken to be the incident laser duration. Now we make some assumptions about the 
maximum value of h. Wilks et al. found that protons gain energy proportional to the electron 
temperature[2], 

Ep = aTe 

where a is somewhere between 2 and 12, depending on the model, and from our simulation results, 
we can roughly estimate a ^ 5. As a result, we can obtain hmax ^ aCgt ~ 2.35(1 + 
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Figure shows the rough range of hole depth, and we can see using our simulation parameters, 
hmin ~ 0.98/im and hmax ~ 4.91/im. It is reasonable that for incident laser with large uq, one 
should use hole-targets with longer ramparts, and the minimum depth shows little change. 


V. CONCLUSION 

In conclusion, proton acceleration in TNSA regime using a tailored hole-target has been inves¬ 
tigated by 2D3V PIC simulations. It is found that proton beams with more intense energy density 
and much smaller divergence angle can be produced from suitably picked hole depth compared 
with those from the single plain target. The dependence of proton beam divergence on depth of 
the hole is also investigated. We give a rough estimation of the depth ranges depending on uq 
of the incident laser, and all the simulation results demonstrate that as long as the depth of the 
hole roughly satisfies the effective action range of the hole target ramparts, the tailored hole-target 
works effectively. 
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